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The Christiansen filter consists of particles of a transparent solid im- 
mersed in a suitable liquid such that the two materials have different dis- 
persions but equal indexes of refraction at some specified wavelength. 
This paper analyzes such a filter as an optical device for providing wave- 
length-dependent resolution. 

We derive the modulation transfer function of this filter by using the 
Huygens-Fresnel diffraction principle. General agreement is obtained with 
experimental measurements made using monochromatic light. 

One proposed single-tube color camera system for picture telephone 
service requires the three primary color component images to have different 
resolutions. An optical filter illustrating the characteristics required by this 
application was designed using the formulas developed in this paper. 
Subjective evaluation of the experimental results obtained with this filter 
indicates no unforeseen degradation of the composite image of the primary 
color components. 

I. INTRODUCTION 

When homogeneous isotropic particles, such as crushed optical glass, 
are immersed in a liquid which has, at a specified wavelength, the same 
index of refraction as the particles but a different dispersion, light at 
the specified wavelength is unaffected in passing through the mixture. 
On the other hand, light of other wavelengths is dispersed. A cell con- 
taining such a mixture of solid particles in an appropriate liquid 
constitutes a Christiansen filter (see Figure l). 1-4 

The wavelength dependent action of the filter can be used to con- 
struct a low-pass spatial filter with a wavelength dependent cut-off. 
In particular, with the proper glass particles and liquid, the filter acts 
as weakly diffusing ground glass for the blue and red images while not 
affecting the green images. 
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An optical device using a color-selective spatial low-pass filter of 
this type is desirable in order to realize a proposed single pickup-tube 
color camera. 5,0 In this camera, color separation is achieved by modu- 
lating the amplitude of the color signals onto different carrier fre- 
quencies. These carrier frequencies are in turn generated by striped 
filters that have different spatial frequencies for the different primary 
colors. It is necessary to reduce the spatial bandwidth of the com- 
posite color image in order to meet the resolution limitation of the 
camera tube. 

This bandwidth reduction can be used because the human eye does 
not readily discriminate between composite images composed of three 
high resolution primary color images, and a composite image com- 
posed of a high resolution green image and low resolution red and blue 
images. 

By considering the filter as a random phase-changing screen, one 
can determine its spatial frequency response as a function of glass 
particle size, thickness of the filter, and position of the filter in rela- 
tion to the image plane. 4,7 

The quantitative understanding thus gained is used to design an 
optical filter with characteristics which are roughly those which 
would be required by an experimental color picture telephone system. 
The experimentally observed effect of this filter indicated no unfore- 
seen degradation of the subjective quality of the final composite image 
transmitted through it. 

n. OPTICAL TRANSFER FUNCTION OF THE FILTER SYSTEM 

Figure 1 shows an optical system transforming an object plane into 
an image plane. If the object plane is incoherently illuminated, the 
system is linear in intensity and therefore: 

i(x, y) = // o(x', y')h{x - x' , y - y') dx' dy' (1) 

where i(x, y) and o(.t', y') are the intensity distribution in the image 
plane and object plane, respectively; and the point spread function h 
of the optical system represents the intensity distribution in the image 
plane for a point object. 
In terms of spatial frequency (f x , /„) response, we have 

IV. , U) = 0(U , U)H(j m , U) (2) 

where the optical transfer function H(f a , f v ) is the Fourier transform 
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Fig. 1 — Optical system with Christiansen filter. 

of the point spread function k. If the transfer function is real and only 
the modulus is used, it is then called the modulation transfer function. 8 
The point spread function /;, by definition, is: 

h(x, y) = u(x, y)-u(x, y)* (3) 

where u(x, y) is the complex amplitude distribution of a point object 
in the image plane. This is called the "point image function," and 
u(x, y)* is its complex conjugate. 

Using the Huygens-Fresnel diffraction principle, the point image 
function u(x, y) can be expressed in terms of the complex field dis- 
tribution A(i, ri) over the exit plane of the filter (see Fig. I). 0-11 We 
should pay particular attention to the fact that we are assuming a 
point object and an isoplanatic region. That is, the optical system is 
assumed linear such that equation (1) holds throughout the entire 
object plane. Thus, for a point object, the point image function is 
given by 

(.r, y) = i~ /£J A& v) oxp j^ [ft - xf + ( V - ?/) 2 ]} r/£ <! v . 

(4) 
The field distribution in the exit plane of a thin Christiansen filter can 
be written as 



MX 



A&n) =^tt^y* u,,> -exp[-^(f + 



•■»] 



(5) 



where $(£, tj) is a random phase function produced by scattering (see 
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the appendix) , and the term E (£, v ) exp[-jk/2d {? + r? 2 ) ] represents 
the spherical wavefront produced by the lens. Since the field outside 
the lens aperture is suppressed, 

wit \ = | constant (£ 2 + 1*)* ^ b 

and b is proportional to the lens aperture as shown in Fig. 1. 

Since the phase function <*>(£, 77) was assumed to be a two-dimen- 
sional random process with isotropic fluctuation, and the experimental 
test object contains only one-dimensional variations in the transmis- 
sion, the problem can be reduced to one dimension for simplicity. Thus, 
by neglecting the constant factor outside the integrals, we may re- 
write equation (4) as 

u(x) = £" BJW™ exp (- I f )-exp {| (« - z) 2 } rff 

-■^jc^^-m* (6) 

The optical transfer function H(j x ) of the system is now obtained 
by taking the Fourier transformation of u{x)u{x)* to give 

/ + CO 
u(x)u(x)* exp (j2irf x x) dx 

= f " B.W* ■ [E& - d\f x ) exp 0*G - d\/,)] j* eft. (7) 
Normalizing (dropping the subscript x) this equation results in 

[* m BJ®t?* w {B& - d\j) exp rj*({ - dXf)]}* (ft 

HO) - ~ = (8) 

/ I E.(Q | 2 (ft 

The statistics of the phase fluctuation *(£) are assumed to be such that 
the average value with respect to £ can be replaced by an ensemble 
average. By denoting an ensemble average by ( ) av , the normalized 
transfer function for this case becomes 

H(j) = DflNexp {*[*(*) - •(« - d\f)}})~ (9) 

where 

Z)(/) = l-f-l. 
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Because the range of the spatial frequency / under discussion is much 
lower than the cutoff frequency produced by diffraction at the lens 
aperture, dkf is much less than 26. Therefore, D{f) can be approxi- 
mated by unity. 

Under the condition that <£(£) is a random gaussian variable of zero 
mean and variance <f> 2 , equation (9) is shown to be 

H(f) =ex P {-$ 2 D [l-p(dX/)]} (10) 

where p{dkf) is the autocorrelation function for the phase fluctuations 
at the exit plane of the filter and is expressed as 

f *({)*({ - d\f) d£ 

P(d\f) = (11) 

j H^ # 

Now, the autocorrelation function P {d\f) must be determined. The 
behavior of the function p is complicated; however, it is known that 

P (d\f) = 1 when dX/ = 
and 

p(d\j) = when d\f is greater than the linear 

dimension of the scattering 
particles. 

In this case there are many possible types of behavior that the func- 
tion P (r) may exhibit. Two common types are. 13-16 

P (r) = exp (-\r/q\) (12a) 

and 

P (r) = exp[-(r/ 9 ) 2 ] (12b) 

where the q represents a correlation length. 

These are both particularly important forms to us because we are 
interested primarily in small r. Thus, equations (12a) and (12b) 
reduce to 



P(r) S* 1 - 
and 



(13a) 



p(r)^l- r - 2 - (13b) 
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The autocorrelation function for any particles, regardless of shape, 
reduces to one of these two forms for small r. In our experiments we 
assumed spherical particles and equation (13a) is found to be the 
more realistic approximation. Thus, let us assume 15,16 



p(d\f) = exp < - 



d\f 



K2a 



(14) 



where the diameter 2a of the particle indicates the order of the size of 
the irregularities. The constant K is a factor expressing the actual cor- 
relation length to the diameter of the particle and can be obtained 
from the function p{dkj) which is empirically determined. Only the 
small values of d\j are of interest so that by approximating equa- 
tion (14) 



p(d\1) S 1 - 



d\f 
K2a 



(15) 



Finally, substituting equations (15) and (27) into equation (10), 
the modulation transfer function obtained is 

J 3.29CL An 2 df \ nc , 

H(1) = exp [- gj j> (16) 

where C is the ratio of particle volume to total volume, L is the thick- 
ness of the filter, and Aw is the difference in refractive index between 
particles and liquid. 

Notice that when / is very large so that p(d\j) = 0, the modula- 
tion transfer function reaches a constant value of 

L/m . **i / 6.5 8CLa An 2 \ , 

H(j) = exp {-*„} = exp«|- - — ^ J- (17) 

This means that there is always finite transmission at high spatial 
frequencies, but that the loss can be made as large as desired by con- 
trolling the parameters of equation (17) . 

III. EXPERIMENTAL 

3.1 Preparation oj Filters 

The filters were made of crown glass type K5 {N d = 1.52320) in a 
solution of ethyl salicylate at room temperature. This combination 
provides a nearly transparent filter in the green region as shown in 
Fig. 2. The dispersion curves were determined by a Carl Zeiss, Model 
A Abbe-refractometer. 
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Fig. 2 — Dispersion curves of optical glass (K5) and ethyl salicylate. 

The glass particles were prepared by pulverizing them with a ball 
mill; they were then cleaned with aqua regia and distilled water. The 
size of the particles was estimated by averaging the dimensions of the 
wire screen used in sorting. The radii of the prepared samples ranged 
from 37/* to 44/*, 44/* to 62/*, 75/* to 125/*, and 125/* to 150/*. The cells 
were made by epoxy-cementing glass plates to 1-inch diameter glass 
tubing of lengths 2mm, 2.5mm, 3.5mm, and 4mm. 

The cells were about half filled with the liquid, then the glass parti- 
cles were poured in slowly so that air bubbles were not earned down. 
The volume of the glass particles was determined from their weight 
and known density of 2.59 grams per cubic centimeter. The volume 
ratio C was found to be approximately 0.42, independent of particle 
size. 



3.2 Measurement of Modulation Transfer Function 

The apparatus used to measure the transfer function of the optical 
system has been described previously. 17 Figure 3 is a rough sketch 
of this apparatus. A target film with sinusoidal intensity distribution 
was used as a test object. The film was imaged through the Chris- 
tiansen filter onto a slit in front of the photomultiplier. Rotation of 
the drum containing the test film caused the image to move across 
the sampling slit. The sinusoidal gratings of the test film varied 
logarithmically from 4 to 300 cycles per inch and were arranged on 
the target film in sequence to allow rapid measurement or direct 
recording of the modulation transfer function. 
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^TRANSPARENCY WHOSE TRANSMISSION VARIES SINUSOIDALLY 
/ AND THE SPATIAL FREQUENCY CHANGE IS LOGARITHMICAL 
ALONG THE DIRECTION OF THE DRUM ROTATION 
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Fig. 3 — Apparatus for measuring the modulation transfer function of the 
filter. 

The measurements were usually carried out in the region of the 
blue mercury arc light to maximize the response of the S-ll photo- 
multipler. The slit width of the photomultiplier was 10G>, which cor- 
responds to a sampling interval equal to one-third the period of high- 
est frequency. An //2.3 lens of 5-cm focal length was used to image 
the target film onto the slit. The target-lens spacing was 6 cm and the 
lens-slit spacing was 30 cm. This resulted in a magnification of 5 and 
gave a range of spatial frequencies from about 0.03 cycles per milli- 
meter to 3 cycles per millimeter. This range of spatial frequencies is 
of interest for potential color picture telephone camera applications. 

The target film was first focused with a cell filled with only the 
liquid placed between the lens and the slit. This compensates for the 
increased optical path length resulting from the nonunity index of 
refraction of the Christiansen filter. 



IV. RESULTS AND DISCUSSION OF THE MODULATION TRANSFER FUNCTION 

The modulation transfer functions of the filters were measured for 
various positions of the filters as well as for various thicknesses and 
particle sizes of the filter. The position of the filter was defined by the 
distance d between the slit plane and the inside of the exit glass plate 
of the cell (see Fig. 1). 

The experimentally determined modulation transfer functions, for 
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several different filter positions with a filter thickness of 2.5 mm and 
particle radius of 100 n, are plotted in Fig. 4. Figure 5 shows the varia- 
tion of the modulation transfer functions with filter thickness for a 
constant particle radius of 100 n and a d of 1/8-inch. From Fig. 4, the 
linear dependence of the exponent 3>„(1 — p) on d can be demonstrated, 
thus validating equation (10). In addition, Fig. 5 shows the approxi- 
mately linear dependence of $ 2 on L, indicating the validity of equation 
(27). From these results, the assumption of equation (15) for the auto- 
correlation function of the phase irregularities also appears reasonable. 
The dependence of the modulation transfer function on particle 
size is shown in Fig. 6. From this figure, the factor K representing the 
ratio of the particle diameter to correlation length [see equation (14) ] 
was obtained and plotted in Fig. 7 as a function of the particle size a. 
A quantitative analysis of the effect of the particle size on the trans- 
fer function is unavailable because the exact expression for the phase 
autocorrelation function in terms of the particle size is not known. It 
is believed that the empirically determined factor K provides a quali- 
tative dependence of the particle size on the phase autocorrelation 
function because only small values of d\f/2aK [see equation (14)] 
are of interest. A theoretical analysis of the phase autocorrelation 
function may not be of great value unless the exact shape of the 
particle is considered. 
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Fig. 4 — Experimental modulation transfer functions for various filter posi- 
tions (o = 100/x, d = 2.5 mm, X = 4358 A). 
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Fig. 5 — Experimental modulation transfer functions for various filter thick- 
ness (a = 100/t, d = 3.18 mm, X = 4358 A). 
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Fig. 6 — Experimental modulation transfer functions for various particle sizes 
of the filter (L = 2.5 mm, d = 3.18 mm, X = 4358 A). 
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Fig. 7 — Factor K of phase correlation length versus particle sizes. 

The linear dependence of the factor K on the particle size a, shown 
in Fig. 7, cannot be extended much further than the radius of about 
150/x because the gaussian phase fluctuation imposes an upper limit 
on a for which the central limit theorem is applicable (see the ap- 
pendix) . In practice, it is not possible to determine the exact correla- 
tion length from Fig. 7, since it is impossible to obtain a large value 
of the correlation length 2aK with sufficient accuracy from extrapola- 
tion of the small values of d\j. 

In order to simulate a situation of independent resolution control of 
primary color images, measurements of the transfer functions at the 
blue, green, and red regions were made with a Christiansen filter 
2.5mm thick and with a 100/x particle radius. In this case, the filter 
was set £-inch away from the image plane. As light sources of primary 
colors, a mercury arc lamp was filtered with Optics Technology, Inc., 
monopass filters, No. 433 for blue (4350 A) , No. 566 for green (5600 
A) , and No. 633 for red (6250 A) .* Figure 8 shows the results. The 
modulation transfer function for the imaging lens itself also is plotted 
in Fig. 8 to show the effect of the lens aberrations. 



♦These narrow band niters were used in order to give quantitative agreement 
with theory. In its proposed application to resolution control the true "wide 
band" primaries should be used. The evaluation of this filter for resolution con- 
trol in a TV color camera is more involved and is the subject of the following 
sections. 
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Fig. 8 — Experimental modulation transfer functions for different wavelengths 
and for imaging lens with a — 100/*, L = 2.5 mm, and d = 3.18 mm. 



Figure 9 shows the comparison between the experimental results 
at different wavelengths and the theoretical modulation transfer func- 
tion given by equation (16). The values of An used in this equation, 
1.5 X 1G" 2 at a = 4350 A, 10" 3 at A = 5600 A and 6 X 10" 3 at A = 6250 
A, were obtained from Fig. 2. The factor K of 1.82 obtained from Fig. 
7 and the measured value C of 0.42 are also used in equation (16) . The 
experimental results are plotted from Fig. 8 after correcting for the 
effects of the aberrations caused by the imaging lens. This correction 
becomes necessary at the higher spatial frequencies encountered with 
red and green light. 

The deviation between the experimental and analytical results 
probably result from thermally induced changes in the refractive 
index of ethyl salicylate, which has a thermal coefficient of index of 
refraction of about 5 x 10-* per C°. A slight change in An causes a 
significant change in the analytical results for the green region, since 
An for green is a small value. On the other hand, the effect of a small 
change in An has very little effect on the calculated blue or red 
responses. Choosing the correct value of An to be used in equation (16) 
is further complicated by the 100 A spectral bandwidth of the primary 
colors used. 

The theoretical and experimental results shown in Fig. 9 agree sub- 
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stantially considering the many approximations made. These ap- 
proximations include the representation of the actual particles by 
spheres with no variation in particle size and assuming a filter thin 
enough to allow the difference in path length through the filter between 
normally and obliquely incident rays to be neglected. 

V. APPLICATION OF THE FILTER 

As Section IV shows, the cutoff frequency of the modulation transfer 
function (spatial frequency response) of the Christiansen filter can 
be controlled by the filter thickness, the particle size, and the posi- 
tion of the filter with respect to the image plane. Using this quantita- 
tive understanding, a filter for independently controlling the resolu- 
tions of the three superposed primary color images was designed. The 
operation of the filter and its suitability for the intended application 
were then experimentally verified. 

In order to evaluate the filter in its proposed application to the 
color picture telephone system, the true NTSC primary color "taking" 
filters were used. 

5.1 Evaluation oj the Filter 

For the color camera, the desirability of having the shortest pos- 
sible distance between the filter and the image plane (that is camera 
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Fig. 9 — Comparison of experimental and theoretical modulation transfer 
functions for selected wavelengths with a = lOO^u, L = 25 mm, and d = 3.18 mm. 
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tube target) dictates choosing a glass particle with about a 40/x radius. 
The filter had a circular cross section 20 mm in diameter, and it was 
2.5 mm thick. 

The experimental arrangement used is shown in Fig. 10. Incandes- 
cent (white) light illuminated a ground glass diffuser which was 
positioned in front of the 35 mm color slide used to simulate a real 
object. In order to be able to separate the three primary color images, 
the incandescent (white) light was actually derived by adding the 
primary colors with the half-silvered mirrors as shown in Fig. 10. 
The primary colors were generated by projecting tungsten filament 
light through the appropriate Wratten filters: No. 25 for red, No. 
47B for blue, and No. 58 for green. 

A Bell and Howell //1.3 Angeniux television camera lens with a 15 
mm focal length was used as an object lens. The image size was about 
10 mm by 10 mm to simulate the raster size of a typical picture tele- 
phone camera. For viewing, the image formed at the image plane 
(that is, vidicon tube target) was magnified about ten times by a sec- 
ond lens. This magnified view of the projected image was used to 
evaluate the filter. 

The Christiansen filter was positioned in front of the image plane; 
the distance between the exit plane of the filter and the image plane 
was adjusted to obtain the desired spatial frequency response. 

In color picture telephone operation, a figure for picture bandwidth 
has not yet been determined. However, at the present time, approxi- 
mations of 1 MHz for the green, 500kHz for the red, and 300kHz 
for the blue images may be assumed for a single pickup camera sys- 
tem. By assuming a horizontal line scan time of 100 f*s and a raster 
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Fig. 10 — Experimental arrangement for evaluation of the Christiansen filter 
as a resolution controlling device. 
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size of 12 mm (horizontal), the relationship between the image res- 
olution and required electrical bandwidth is determined. 

The spatial frequency response was measured by using a pie-chart 
test target disk with 30-cycle gratings along the circumference. Since 
the spatial frequency of the grating varies continuously along the 
radial direction of the disk, the cutoff frequency of the image can be 
measured by knowing the point where resolution fails. The image 
size of the test target was 7.5 mm in diameter. 

In order to determine the effectiveness of the filter as a resolution 
controlling device, an image without the filter was viewed for com- 
parison. In this case the Christiansen filter was replaced by a cell the 
same size as the filter, filled with only the ethyl salicylate solution. 
This maintains the same optical path lengths within the imaging 
system. 

5.2 Image Quality through the Filter 

In order to determine the resolution, each pictoral subject imaged 
through the system was associated with the corresponding image of 
the resolution test pattern. This was done by replacing the object slide 
with the test target object when photographing the system output. 
The resolutions were determined simply by observing the images of 
the test pattern resolution chart. The spatial frequency of the test 
image was 1.27 cycles per mm at the outer edge of the circular pat- 
tern. This spatial frequency corresponded to 153 kHz by assuming 
100 microseconds for the line scan time and 12 mm for the raster size 
of the camera tube. 

Figure 11a is a photograph of the image formed by incandescent 
(white) light without the Christiansen filter; Figure lib is the cor- 
responding photograph of the image of the resolution test pattern. 
Notice that the same quality images were obtained with separate 
illumination of each primary color since achromatic lenses were used 
in the experiment. 

Figure 12a is a photograph of the image formed by incandescent 
(white) light with the Christiansen filter placed 1 mm in front of the 
image plane. Figure 12b is the corresponding image of the resolution 
test chart. The distance of 1 mm between the filter and the image 
plane was chosen to demonstrate the particular case for which the 
filter will be used in the proposed color camera system. Some of the 
dark spots in the photographs were caused by the presence of dirt 
in the filter. Blurring of the image by low-pass filtering is noticeable. 
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Fig. 11 — (a) Portrait and (b) resolution test chart showing images formed by 
incandescent (white) light without the Christiansen filter. 

This was introduced deliberately by the Christiansen filter and should 
not be considered as a defect of the optical system. The quality of the 
image seems acceptable in all other respects. 

In order to illustrate to what extent each primary color image in 
the composite image (Fig. 12) is spatially filtered, three separate 




Fig. 12 — (a) Portrait and (b) resolution test chart showing images formed by 
incandescent (white) light with the Christiansen filter inserted in the optical 
system. 
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photographs were taken with each primary color by blocking the other 
two color sources. Figures 13a and b show photographs of the images 
formed with green primary light only. The image would be as sharp 
as that in Fig. 11 if the filter was completely transparent at green. 
However, the spectral bandwidth of the NTSC green primary color is 
wider than the matched bandwidth of the filter in green, so that the 
quality of Fig. 13 is not quite that of Fig. 11. This illustrates that it 
is advantageous to limit that bandwidth of the green image so that 
it does not extend beyond the required minimum value. 

Figures 14a and b show the images formed with red primary light 
only. Blurring of the image caused by low-pass spatial filtering by the 
Christiansen filter is noticeable. The cutoff frequency of the image, 
determined from Fig. 14b, is about 4 cycles per mm (the resolution 
falls about one-third of the radius of the pattern) corresponding to a 
500 kHz electrical signal. 

Similarly, Figs. 15a and b show photographs of the images formed 
with blue primary light. Blurring of the image is much more notice- 
able than that of Fig. 14, as expected from characteristics of the 
Christiansen filter. The cutoff frequency of the image is about 2.5 
cycles per mm (the resolution falls about half of the radius of the 
pattern) corresponding to a 300 kHz electrical signal. 

It should be emphasized that although the blurring of the red and 
blue images is very noticeable when they are viewed alone (as in 
Figs. 14 and 15), the composite image of all three primaries is of 
much better quality. This indicates that a significant amount of the 
detail in the composite image is provided by the green component. 18 




Pig. 13 — Images formed by green primary component color with the filter 
inserted in the optical system. 
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Fig. 14 — Images formed by red primary component color with the filter 
inserted in the optical system. 

In order to illustrate the subjective quality of a composite color 
image transmitted through this system, Figs. 11a and 12a are printed 
in color as Figs. 16 and 17, respectively. These photographs were taken 
on Ektacolor L film. 

Much work remains to be done to determine precisely what resolu- 
tion values should be used, but such a study is beyond the scope of 
this paper. The point we are trying to make is that we can control the 
resolution of the three primary color component images and that there 
is no appreciable degradation beyond that deliberately introduced in 
the form of spatial bandwidth reduction. Although this system does not 
have a sharp cutoff characteristic, that is, it does not have a rectan- 




Fig. 15 — Images formed by blue primary component color with the filter 
inserted in the optical system. 
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gular passband fsee equation (16)], it can nevertheless reduce the 
image detail of a selected color component by any desired amount, 
while still permitting coarser image structure to remain. 

VI. CONCLUSIONS 

Formulas for the modulation transfer function of the Christiansen 
filter have been derived. The filter thickness, particle size, position of 
the filter from the image plane, and difference between the particle 
and liquid indices of refraction have all been shown to control the 
cutoff frequency but not the shape of the modulation transfer func- 
tion. The modulation transfer functions as measured with essentially 
monochromatic sources were in quantitative agreement with those 
predicted from theory. 

The proposed use of the Christiansen filter to reduce the red and 
blue color detail in a composite image is illustrated using the true 
"wide band" primary colors. Subjective evaluation of the composite 
image shows that there is no appreciable degradation beyond that 
deliberately introduced in the form of spatial bandwidth reduction. 
The result confirms the applicability of the filter in color picture tele- 
phone-like schemes as a means of independently controlling the spatial 
resolution of the component color images. 
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APPENDIX 

Random Phase Variation in the Chmtia?isen Filter 

When a plane wave is normally incident on the entrance plane of 
the Christiansen filter, the wave emerges from the filter with the phase 
varying across its wave front but with no appreciable change in ampli- 
tude. The particle size is on the order of hundreds of a wavelength; 
the difference in refractive index between the particles and liquid is 
a few parts per thousand. Therefore, the distribution of phase over 
the emerging wave-front can be described approximately in terms of 
the distribution of the path length of each ray through the particles. 

Since the particles are close-packed, it is assumed that the average 
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number of particles traversed by a ray as it crosses the filter is uni- 
form. The particles are considered to be spherical with constant radii. 
The statistical fluctuation in the phase across the exit plane of the 
filter is caused by the ray passing through different parts of the 
particles. 

Assume the ratio of particle volume to total volume is C and the 
cross-sectional area of the filter is S. The average number of particles 
intersected by a single ray which passes through a filter of the thick- 
ness L can then be calculated as 



N = 



__ ( C8L \(^\ - ? QL 



W A s 



(18) 



where a is radius of the particle. The first parenthesis represents the 
total number of particles in the filter. 

If the ray is traveling in the z direction in the Cartesian coordinates 
system (x, y, z), the ray path length through an individual particle 




Fig. 16 — Ektacolor print of Fig. 11a. 
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Fig. 17 — Ektacolor print of Fig. 12a. 

I will be a function of the lateral displacement of the ray from the 
particle center in direction .r and y. The average path length through 
a particle is 

(l)„ = f f" ' 2(a 2 - x 2 - if) h P(x, y) dx dy (19) 

where P{x, y) is the probability density function. The lateral dis- 
placements in directions x and y are uniformly distributed random 
variables in the range of the area war. Therefore, P(x, y) is l/va* and 
equation (19) becomes 

</)av = f r 2(« 2 " x * ~ y 2 ) 1 -^dxdy = U- (20) 

This result differs from that obtained by R. H. Clarke, whose prob- 
ability density function (equation A4 of Ref. 4) depends upon one of 
the spatial coordinates. 
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The second moment of I is 

(l\ v = f" f 4(a 2 - x 2 - if) — 2 dx dy = 2a 2 . (21) 

Therefore, the variance of the random variable / is 

<Oav - <0 2 av = 2a 2 - (ia) 2 = |a 2 . (22) 

If the plane wave E„(x, y) (traveling in direction z) is normally 
incident to the entrance plane of the filter, the field distribution at 
the exit plane of the filter can be expressed by the equation 

E(x, y) = E (x, y) exp {;*(*, y)\ (23) 

where #(.r, y) is a random phase function. 

The phase of the rays at the exit plane of the filter with reference 
to the average phase across the plane can be written as 

* = E k An(l m - (Oav) (24) 

m = l 

where k = 2-n-Jk (a is a wavelength in free space) and An is the differ- 
ence in refractive index between particles and liquid. The process 
corresponding to each m is an independent random variable so that 
the summation according to the central limit theorem should be gaus- 
sian. This gives 

p($) = (2^: exp( "* 2/24>2o) - (25) 

Thus, the variance of the phase at the exit plane is 

*l = ([S fc An(/. - (U«)J) = N(k An) 2 ((l - (Z) av ) 2 > av 

+ N(N - l)(fc An) 2 <(/ - (Z> av )> 2 av . (26) 

Noting that the second term is zero and using equations (18) and (22), 
equation (26) becomes 

*2 XT/? a \2 2 2 6.58CLa An /0 _ N 

3>o = N(k Any fa = -^ (27) 
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